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Abstract 

It is demonstrated that the complex scahng method can be used in practical 
calculations to localize three-body resonances. Our model example empha- 
sizes the fact that in three-body systems several essentially different asymp- 
totic behaviors can appear. We show that the possibility of these different 
asymptotic configurations can lead to an apparent, resonance like structure 
in the three-body continuum. 
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Three-body resonances, systems which decay into three-body final states, play important 
role in few-body physics. The description of these states are somewhat easier than that of 
the scattering processes which lead three (or more) particles in the outgoing channel. The 
study of three-body resonances hopefully will help us to clarify some points of the still not 
totally solved scattering problem. 

One of the main difficulties in the general formulation of the quantum mechanical N-body 
resonance (and scattering) problem is the prescription of the asymptotic behavior. There is 
fundamental difference between the problems where only short range interactions occur and 
those where a long range force, e.g. the Coulomb force, is present [Q. There are methods 
which eliminate the explicit reference to the unknown or partially known asymptotics, e.g. 
the J-matrix method 0, and the potential separable expansion method P], but up till 
now they have been developed only for two-body systems. In the case of bound states the 
question of the asymptotic behavior is not a serious problem. The majority of the methods 
work because in most cases the bound state asymptotics hardly affects the determination of 
physically observable quantities. However, in the case of resonances and scattering states, 
the asymptotic behavior plays a crucial role. 

The complex scaling method (CSM) ||^ reduces the description of resonant states to 
the description of bound states, thus avoids the problem of asymptotics. This method 
handles the non-Coulomb and the Coulomb cases on equal footing. There is an extension 
of this method to two-body scattering states, too 0. It is an intriguing question whether 
this method can be extended to the general many-body scattering states. Up till now, 
three-body resonances have been investigated e.g. using the Faddeev method 0, by real 
stabilization |]^, and in the time delay matrix formalism In this paper we study them 
using the complex scaling method. Although the CSM was used to describe three-body 
systems above the three-body breakup [0, the authors did not specify how they identified 
three-body resonances. This research was prompted by the contradicting theoretical results 
concerning the existence of a soft dipole resonance in the neutron halo nucleus ^He []TD|,|TT|. 
What makes this problem difficult, is the mixture of the underlying nuclear physics and 
three-body dynamics. Here we want to clarify some points of the three-body dynamics. 

For the sake of simplicity, we recall the main points of the CSM in two-body context. 
In the coordinate space, resonance eigenf unctions, corresponding to the complex energy 
solutions of the 

H\^iJ) = {f + V)\^) = E\^) (1) 

Schrodinger equation, show oscillatory behavior in the asymptotic region, with exponentially 
growing amplitude, ~ exp[i{K — i'y)r] (k, 7 > 0). Thus, they are not elements of the space. 
The complex scaling method means that instead of Eq. (|1|), we solve the eigenvalue problem 
of the transformed Hamiltonian Hg = U{6)HU^^{6): 

Hel^e) = Eel^e). (2) 

U{6) is an unbounded similarity transformation [jl2|, which, in the coordinate space, acts 
on a function /(r) such that 

U{9)fir)=e'^'/'f{re^'). (3) 
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(If 9 is real, U{9) means a rotation into the complex coordinate plane, if it is complex, it 
means a rotation and scaling.) The two problems are connected by the Aguilar-Balslev- 
Combes theorem If \^ is a (dilation) analytic operator, then (i) the bound eigenstates 
of H are the eigenstates of Hq, regardless of the actual value of 6*, within < 6' < 7r/2; (ii) 
the continuous spectrum of H will be rotated by an angle 26*; (iii) a complex generalized 
eigenvalue of Eq. (H), i?res = s — i\T, e, F > (with the wave number fcres = n — ij, 7 > 0), 
belongs to the proper spectrum of He provided 29 > |argi?i.es|- Roughly speaking, the 
complex scaling transformation changes the asymptotic wave function from exp[i{K — i'j)r] 
to exp[i{K — i'y)r exp{i9)], which, in the case of 29 > |argi?res| = 2| arg fcresl, makes the 
diverging wave function localized. 

If we have particles, we can transform the problem from one-particle coordinates to 
certain interparticle relative coordinates (Jacobi coordinates): 

{ri, r2, . . . , Tn} {ti, ta, . . . , tAr_i}, (4) 

where the origin is fixed at the center of mass, so that J2i^=i "^i^^i = 0, where mj are the 
particle masses. The application of the CSM in this case means that we transform all relative 
coordinates under the action of tJ{9). Note that the (^) transformation is linear, therefore 
if we carry out the complex scaling transformation in a certain set of Jacobi coordinates, it 
results in the same transformation in all other possible sets of relative coordinates. It means 



that we cannot choose the rotation angles independently, like in the multichannel CSM [|14 
in the different configurations. 

Our present model problem consists of three particles with masses mi,m2, and 7713. 
We choose mi = 1112 so as to have only two different systems of Jacobi coordinates. The 
interactions between the particles are separable forces 

% = \^Q{h))\j{ipo{h)l J > 1 = 1,2,?, (5) 

where |<^o(^)) is the eigenfunction of the three dimensional harmonic oscillator with n = I = 
0, b is the oscillator size parameter, and Xij are the potential strengths. In coordinate space 
the interactions depend on the various relative coordinates. Each interaction has a natural 
Jacobi coordinate system in which it has the simplest form, e.g. in the (12)3 coordinate 
system (there is a relative coordinate between particles 1 and 2, ti2 and another between 
(1,2) and 3, t(i2)3) the Vu interaction depends only on ti2. 

The kinetic energy operator is easy to express in any coordinate system, e.g. in (12)3 it 
looks like 



T = 

2 



J_ 1 

^ti2 + ^t(12)3 
>il2 /i(12)3 



(6) 



where the Laplace operators are differential operators in the appropriate Jacobians. It is 
easy to show that the application of the (^ complex scaling transformation to our potential 
and kinetic energy operators is equivalent with the change of b to bexp{i9) in (^), and the 
multiphcation of the right-hand side of (|^) by exp{—2i9), respectively. 

As the CSM localizes the resonant wave function, we can use any bound-state method 
to describe them. Here we use the wave function expansion method. We consider three 
different wave functions. In coordinate space 
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^1 = Z!c»i<^i(^12)¥'i(t(12)3), (7) 
^2 = 5Z'^uV'i(*23)V'i(t(23)l), (8) 

and the unification of tlie above two, \E'3 = + \I/2. Tlie expansion coefficients c and d 
are to be determined from a variational principle. Our potentials act only between s-waves, 
therefore each oscillator function carries zero angular momentum in (|^) and (P). We choose 
the oscillator size parameter of the wave functions, 6, different from h in order to make the 
trial functions more flexible. The summation limits in the wave functions are chosen to 
reach stable convergence. 

For the necessary matrix elements we need to calculate the overlap of the product oscil- 
lator states between different Jacobi coordinate systems, and the Laplace operator between 
such states. As examples, we show typical terms that appear in the overlap, kinetic energy, 
and potential matrix elements: 

(V^i(t23)V5i(t(23)l)|V5fc(tl2)V5z(t(12)3)), (9) 



(<^i(t23)<^j(t(23)l)|At,J(pfc(ti2)v3;(t(12)3)), (lO) 

and 

(V9i(t23)v2i(t(23)l)|V5o(tl2))(V5o(tl2)|V5fe(tl2)V5Kt(12)3))- (H) 

Using the Talmi-Moshinsky-Tobocman transformation [|15| we can express a product of os- 
cillator states, given in a certain Jacobi coordinate system (a), in terms of product oscillator 
states in another system (a'), e.g. 

V,nt23)¥^- (t(23)l) =T.<Ukit,2M\Hl2)z). (12) 



where the sum is finite, and the transformation coefficients can be calculated, e.g. by the |T6 
program. Using these transformations, and in addition the overlap between two oscillator 
functions with different size parameters, and the matrix element of the Laplace operators be- 
tween such oscillator functions (for the formulae see e.g. Jl^); all necessary matrix elements 
can be calculated analytically. 

We choose mi = m2 = 2, and = 4 = 1 and atomic mass units are used), 
h = 1.0, and b = 2.0. The use of separable interactions allows us to set up their strengths 
in such a way that resonances occur at prescribed energies in the two-body subsystems 
T^ . The choice A12 = 0.6377 + ^0.0697 results in a resonance in the (1,2) subsystem at 



E = 1.5 — iO.5. The A13 = A23 = 1.0 strengths give a resonance in the (1,3) and (2,3) 
subsystems at 1.7553 — ^0.2438 energy. As an illustrative example we show in Fig. |l|(a) the 
result of a GSM calculation for the (1, 2) subsystem. The working mechanism of the method 
is clearly seen, the discretized continuum points are rotated, and the resonance is revealed. 

In Figs. |I](b)-(d) we show the results of the three-body GSM calculations using the 
\l/2, and \l/3 trial functions, respectively. We can see that as the rotation angles are large 
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enough to localize the resonances in the subsystems, there are discretized continuum points 
lying on strait half-lines which start from the position of the resonances of the subsystems. 
These starting points act as non real thresholds. For example, in Fig. |l](b) the half-line starts 
at 1.5 — iO.5 which is the resonance energy in the (12) subsystem. This is in full agreement 
with the mathematical theorems p9|j20|] . In addition to these continuum points, we can see 
that there is an isolated point at 4.128 — 20.337 in each figure. We can identify this point 
as a three-body resonance. The fact that this point occur in each figure shows that this 
state of the three-body system can show up both 3(12) and 1(23) asymptotic behaviors. 
This is exactly what we expect from a three-body resonance. This behavior gives us a 
method to identify three-body resonances in practical calculations. But what can we say 
about the continuum states which lie on the half-lines starting from the resonance energies 
of the subsystems? In these states e.g. in the 1(23) configuration there is a resonant state 
in the (2,3) system and a scattering state between 1 and (2,3). These continuum states 
are essentially different from a pure three-body scattering state which is represented by the 
continuum points lying on the half-lines start at the origin. The resonance+scattering type 
continuum states can represent some kind of sequential decay, where the life time of a quasi- 
stationary subsystem, (2,3), is longer than the time needed for 1 and (2,3) to be scattered 
off. 

The fact that different type of continuum states can be present in a three-body system 
means that the three-body continuum has structure in addition to the three-body resonances. 
Let us speculate a bit about this additional structure. In Fig. |] we show the distribution 
of the continuum points (in 0.2 wide energy bins) of the model whose wave function is \l/3. 
We note here that there is no qualitative difference between not complex scaled {6 = 0) 
and complex scaled results because the CSM causes only a contraction of the continuum 
points, while the resonances remain stable. We can see in Fig. ^ the resonant structure 
around 4.2. But, in addition, we see structures, the concentration of the continuum points, 
around 1.6 and 2.0 which values coincide the real parts of the resonance energies of the 
subsystems. These resonant-like structures are apparent and they are the consequence of 
the fact that if the energy is larger than the threshold energy of a subsystem's resonance, a 
new, resonance+scattering, asymptotic behavior can appear. We should note that the shape 
of the background distribution in Fig. ^ is surprising. In the case of two particles, the wave 
function expansion method can be considered as if we closed our system into a box whose 
size is the finite spatial region of the trial function. This leads to a ~ {n = 1,2,.. .) 
spectrum. From this, it follows that the number of continuum points which are in a AE 
interval around E (if AE is small) is AA^ ~ 1/ \fEAE. We checked in a two-body model that 
this is a really good approximation. Our three-body spectrum in Fig. 2 is, however, strongly 
differs from such a shape. In the case of three particles, the Schrodinger equation can be 
cast, in the hyperspherical coordinate, into a form which is similar to a two-body equation, 
see e.g. Ref. ||21[]. If all angular momenta are zero, in this reformulated Schrodinger equation 
a centrifugal barrier occures with L = 3/2. This nonzero L really leads to a reduction of the 
low energy solutions, but this is a marginal decrease. There must be another effect which 
supresses the low energy spectrum. 

In conclusion, we showed that the complex scaling method can be used in practical 
three-body calculations. In this method, three-body resonances can be identified by those 
resonant energy solutions which appear in all Jacobi coordinate system. We pointed out that 
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the possibility of resonant+scattering type asymptotic behaviors can lead to an apparent 
structure in the three-body continuum. The case of the ^He soft dipole mode is very similar 



to our present example. That nucleus is a genuine three-body, a + n + n, system p2| , p3 
There are two resonances in the a+n subsystem at 0.89 — iO. 60 MeV and 4 — i4 MeV energies 
p^ . (And there is an antibound state in the n+n subsystem, which is out of the scope of this 
article.) From what we can learn our present model, it is fairly possible that the structure 
in the ^He continuum, which is interpreted as a fingerprint of a three-body resonance, the 
so-called soft dipole resonance, is nothing but a consequence of the three-body dynamics. 
Of course here we considered only one side of the problem. The question of the ^He soft 
dipole resonance can only be answered in a model which contains both the nuclear physics 
and the three-body dynamics properly. The application of the complex scaling method in a 
realistic model of ^He is in progress. 
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FIGURES 



FIG. 1. Energy eigenvalues of a complex scaled (a) two-body problem with mi = m2 = 2, 
and Ai2 = 0.6377 + i0.0697; (b)-(d) three-body problem with rrii = m2 = 2, = 4, 
Ai2 = 0.6377 -|- i0.0697, and A13 = A23 = 1.0. The trial wave function is: (b) ^'i, (c) ^'2, and 
(d) ^'3. The 9 rotation angle is 0.4 rad in each figure. 

FIG. 2. Distribution of the continuum energy solutions of the three-body problem specified in 
Fig. |l], with the trial function ^'3. The solutions are groupped within 0.2 wide energy bins. 
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